Characterizing quantum phase transitions by single qubit operations 
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We introduce observable quantities, borrowing from concepts of quantum information theory, for 
the characterization of quantum phase transitions in spin systems. These observables are uniquely 
defined in terms of single spin unitary operations. We define the energy gap between the ground state 
and the state produced by the action of a single-qubit local gate. We show that this static quantity 
involves only single-site expectations and two-point correlation functions on the ground state. We 
then discuss a dynamical local observable defined as the acceleration of quantum state evolution 
after performing an instaneous single-qubit perturbation on the ground state. This quantity involves 
three-point correlations as well. We show that both the static and the dynamical observables detect 
and characterize completely quantum critical points in a class of spin systems. 

PACS numbers: 03.67.-a, 75.10.Jm, 05.50.+q 



Quantum phase transitions differ from classical ones 
in that they are driven by quantum fluctuations. They 
are induced by the changes in external control parame- 
ters or coupling constants, which can modify the ground 
state of a quantum system [l| . In analogy with classical 
phase transitions driven by thermal fluctuations, quan- 
tum phase transitions are usually analyzed in terms of 
the scaling behavior of equilibrium properties and diverg- 
ing correlation lengths as signatures of quantum critical 
points X^J- I^ ^1^*5 l^^t years interest in quantum critical 
phenomena has increased enormously thanks to several 
studies that have highlighted some intriguing connections 
between quantum phase transitions and the scaling be- 
havior of nonlocal quantum correlations in the ground 
state of complex quantum systems |3|, jj, |5|. Quantum 
critical points in some classes of quantum spin models can 
be characterized by the divergence in the first derivative 
of the two-site (two-qubit) concurrence, a proper measure 
of bipartite entanglement [a| , with respect to the proper 
control parameter, even though recent work on different 
systems has revealed that this behavior of the concur- 
rence is not universal at generic quantum critical points 
[21 la. Is! I ^^'^ different classification schemes have been 
proposed recently jlfll, |ll|. The studies on the scaling 
behavior of entanglement prompt naturally the question 
on the possible links of entanglement measures to observ- 
able quantities in quantum critical systems. More gener- 
ally, as quantum entanglement plays a fundamental role 
both in quantum information and quantum computation 
|l2l |. the further question arises whether it is possible to 
reveal structural aspects of quantum critical systems us- 
ing other concepts and techniques, beyond entanglement, 
motivated and/or imported from the arena of quantum 
information theory. 

Following this line of thought, in this work we discuss 
the behavior of observables associated to unitary trans- 
formations on individual spins, i.e. local unitary single 
qubit operations. We will show that they can be used to 



detect and characterize quantum critical points in a class 
of quantum spin systems. First, we will focus our atten- 
tion on the energy gap needed to realize single qubit oper- 
ations by local excitations above the ground state. Next, 
we will consider the short-time diffusion coefficient of the 
on site magnetization along a given direction, i.e. the 
quantum speed characterizing the dynamics of the mag- 
netization under the action of single qubit operations. 
We will first establish that the observable quantities as- 
sociated to fundamental local operations (elementary lo- 
cal gates) can be unambiguously expressed as functions 
of one-, two-, and three-point correlation functions eval- 
uated in the ground state. These relations in turn will 
allow to study their behavior when approaching a quan- 
tum phase transition, and to detect and characterize it. 
Let us consider the class of spin- 1/2 ferromagnetic sys- 
tems described by the Hamiltonian 



H 



4e[(i 



7KV+i + (l-7KV+i]+/iE^^' 



(1) 

where i runs over the sites of the lattice, J is the nearest- 
neighbor coupling constant, h is the transverse magnetic 
field directed along the z-axis, 7 is the anisotropy pa- 
rameter, < 7 < 1, erf are the Pauli matrices (a — 
X, y, z), and periodic boundary conditions are assumed 
throughout. It is convenient to define the dimensionless 
reduced couplin g co nstant A = J /h. These models are 
exactly solvable |l3lll4Lll5|. and, in principle, the struc- 
ture of the correlation functions in the ground state can 
be studied accurately. For 7=1 Eq. Q reduces to the 
Ising model whereas for 7 = it reduces to the isotropic 
XY model. In the interval < 7 < 1 the models belong 
to the Ising universality class. For < 7 < 1, in the ther- 
modynamic limit the system undergoes an order-disorder 
phase transition at the critical value Ac = 1. If A > Ac 
the magnetization < a'f > is finite and vanishes at the 
transition, while the magnetization < af > remains finite 
for any value of A and vanishes only in the limit A — > 00. 



The case 7 = is in a critical regime for all A > Ac- 

A local gate (or one qubit gate) is any unitary transfor- 
mation acting as the identity 1 on any subspace defined 
on any spin of the lattice except one. Consider the single 
spin, on which the local gates act nontrivially, located, 
say, at a certain site k. Then, the corresponding unitary 
operator can be written as U ^ Yii^k ^iOk, where Ok 
is a generic_unitary operation acting on the spin placed 
on site k 12] . The natural basis in the space of local 
gates at lattice site k is formed by the three Pauli ma- 
trices cr^, tJ^, (7|, and the identity 1^. For instance, the 
Hadamard gate is if^ = [uf, + (t|)/\/2, the phase gate is 
Sk ~ ((1 — i)<^k + (1 + *)lfc)/2, and the spin flip operator 
is simply cr^. 

The action of a local gate on the ground state \g) 
of the system takes it to a state |V') — Ok\g) that, in 
general, is not an eigenstate of the many-body Hamil- 
tonian H. Introducing the (dimensionless) energy gap 
AE{Ok) = HMHW - {g\H\g))/h it is possible to de- 
termine the exact expressions for the three energy gaps 
associated to the fundamental gates in terms of various 
correlation functions. We obtain: 



AS, = 2A(1 - -i)Gyy - 2M, , 
ASj, = 2A(1+7)C. -2M, , 
AS, = 2A [(1 + 7)C + (1 - -f)Gyy\ , 



(2) 



where ASq, denotes the energy gap due to the action of 
the gate tr^, M, is the value of the on site magnetization 
along the z-axis in the ground state: M, = {g\af.\g), 
and Gaa denote the values of the nearest neighbor two- 
point correlation functions in the ground state: Gaa ~ 
(5|cr^cr?j,_|_^Jg). In the derivation of Eqs. Q we properly 
exploited the invariance of Hamiltonian Eq. |^ under 
translations. 

The major significance of the set of Eqs. (0) is that 
the energy gaps between excited states obtained by local 
gate operations on the ground state and the ground state 
itself can be expressed as functions of (bipartite) expec- 
tation values on \g). Therefore, correlation functions and 
other ground state properties can be determined once the 
energy gaps entering in Eqs. ^ are known, and changes 
in \g) must be reflected in changes of the energy gaps. 

In FigQ] we plot the behavior of the derivatives with 
respect to A of the three energy gaps in the case of the 
Ising model (7 = 1). Here and in the following of this 
work we will exploit the studies performed by Barouch 
and McCoy [lj| to obtain exact analytical expressions. 
We find that the derivatives of the three energy gap di- 
verge logarithmically at the critical point Ac = 1 as 



dAEa 
dX 



Cai"t) In |A — Ac| + const 



(3) 



where, for the Ising model, the values of the coefficients 
Ca{l) are c,(l) = -1.276, c,(l) = 0.638, and Cy{l) ^ 
-c.(l). 




1.3 



FIG. 1: First derivatives of the energy gaps as functions of A 
in the Ising model (7 = 1). Red line is dAE^/dX, blue line 
dAEy/dX, and black line ~-dAEx/d\. 



For < 7 < 1 the models described in Eq. Q are 
in the same universality class of the Ising model. With 
respect to the concept of universality - the critical prop- 
erties depend only on the dimensionality of the system 
and on the broken symmetry in the ordered phase - it 
is important to verify that the energy gaps show simi- 
lar behaviors when the models in this range of values of 
7 approach the quantum critical point. We have veri- 
fied the same property of logarithmic divergence when 
the anysotropy parameter varies in the open range (0, 1). 
For instance, for 7 = 0.5 the derivatives of the energy 
gaps close to the quantum critical point obey the same 
law Eq. Q with c^(0.5) = -2.550, Ca;(0.5) = 0.638 and 
Cj,(0.5) — — Ca;(0.5). In general, for decreasing 7, we find 
that the value of c^ decreases and, remarkably, that the 
values of Cx and Cy, and the relation between them, re- 
main constant. 

Considering the isotropic XY model (7 = 0), the be- 
havior of the energy gaps close to the quantum critical 
point changes radically because the system belongs to a 
universality class characterized by a new symmetry, the 
invariance under rotations around the z-axis. This sym- 
metry imposes the plane degeneracy E^ = Ey in Eqs. |(5J. 
In Fig[21 we show the behavior of the derivatives of the 
two remaining independent energy gaps {E^ and E^) as 
a function of A approaching the critical value A^,. One 
finds that the derivative of the energy gap AS, asso- 
ciated to the TT phase gate operation a^ diverges as an 
inverse power approaching Ac from above: 



dAE^ 
dX 



= °(A^-A?)-V2 



(4) 



The derivative of the energy gap AE^ associated to the 
spin flip gate remains finite for all values of A. Hence, 
only the energy gap associated to the (tt) phase gate op- 
eration detects the onset of the critical regime in the XY 
model, due to the larger symmetry constraints. 

Next, we consider quantum state dynamics rather than 
energy observables. In general, a local gate operation 
perturbs the ground state \g), and under the subsequent 
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FIG. 2: First derivatives of the two independent energy gaps 
as a function of A for the XY model (7 — 0). Red line: 
dAE^/dX; black line: dAE^/dX = dAEy/dX. 



dynamics driven by Eq. |^ any observable that does not 
commute with the system Hamiltonian acquires a non- 
trivial time evolution. It is convenient to focus on the 
dynamics of the magnetization along the z-axis on site 
k, because it is nonvanishing for all values of A, and hence 
its time evolution is best suited for analytical treatment 
as well as for experimental observation. Let us then de- 
fine the dimensionless time r — {ht)/h, and consider the 
short-time dynamics (Ar <C 1) of the magnetization on 
site k after the local gate operation has been performed 
at T = 0. One has < o-|(t) > - < cr|(0) >= -AAaiXr)'^, 
where < cr|(0) > is the magnetization along the z-axis 
evaluated on the perturbed state at time 0, < <J^{t) > 
is the on site magnetization at short times t after the 
action of the local gate, and A^ is the short-time diffu- 
sion coefficient, or acceleration, whose form depends on 
the single qubit operation that has been performed. This 
acceleration can be expressed, similarly to the case of the 
energy gaps previously studied, as a linear combination 
of correlation functions evaluated in the ground state 



FIG. 3: First derivatives of the short-time diffusion coeffi- 
cients Ax (black line), —Ay (blue line), and A^ (red line), as 
functions of A for the Ising model (7 = 1). Notice that A^ is 
always constant as a function of A. 



n > 2 is related to the dynamical nature of the observ- 
able. Including next-to-leading order powers of At in the 
expansion of < (T^(r) > — < cr^(O) >, one would need to 
evaluate four-point and higher-order correlations. 

In Fig|21we plot the behavior of the derivatives of the 
diffusion coefficients dA^/dX (a = x, y, z) as functions of 
A for the Ising model (7 = 1). The derivatives of A^, and 
Ay diverge logarithmically approaching Ac, in analogy 
with the derivatives of the energy gaps (see Eq. Q and 
FigP), as 



dAg 
dX 



= da (7) In |A — Ac| + const 



(6) 



The coefficients da (7) can be evaluated analytically. In 
the Ising case the following symmetry holds: dx{l) — 
— Ca;(l), dy{l) = —Cy{l). lu thc wholc interval < 7 < 1 
the accelerations A^^ and Aj, exhibit a very similar scal- 
ing behavior in the proximity of the quantum critical 
point, obeying in all cases Eq. ©. The coefficients 



For the three fundamental operations (ct^, ct^, and cr^^.), d^Al) in Eq. © are related to the coefficients 0^,^(7) in 



the coefficients A^ can be worked out analytically and 
read, respectively. 
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Comparing Eqs. (|2Jl and Eqs. (jSJ, we see that in the lat- 
ter there appear three-point correlations involving near- 
est neighbor and next to nearest neighbor sites: Gapa = 
{9\'^k-i'^k'^k+i\fi) ■ ^^^ presence of correlations of order 



Eq. |(31l by a simple translation of minus the latter by an 
amount that decreases monotonically with the value of 
the anisotropy and vanishes when the latter reaches its 
maximum 7=1. As an example, at 7 = 0.5 the coeffi- 
cients 4(0.5) and dy(0.5) are 4(0.5) = -c^^ (0.5) -1-0.239 
and dyiQ.h) = -Cy{Q.h) -f 0.239. 

The reason why dA^/dX does not exhibit scaling in the 
Ising model can be understood by comparing with its be- 
havior at different values of the anisotropy. One finds 
that dAz/dX scales logarithmically, just hke dA^/dX and 
dAy/dX, when approaching the quantum critical point in 
the whole interval < 7 < 1. At first sight, this differ- 
ence in the behavior of the acceleration A^ for the Ising 
model versus the other anisotropic cases would seem to 
imply a violation of the universality principle. However, 
studying the behavior of A^ as a function of A when the 
anisotropy parameter 7 approaches the Ising limit, one 
finds that dA^/dX still obeys Eq. ©, but with a coef- 
ficient (^2(7) that decreases monotonically and vanishes 
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FIG. 4: Behavior of the first derivatives of the velocities A^ = 
Ay and of the velocity Az as functions of A when approaching 
the critical point Ac for the isotropic XY model (7 = 0). Red 
line: dA^/dX; black line: dA^/dX = dAy/dX. 



exactly at 7 = 1. 

In the case of the isotropic XY model, we enter in 
a different class of universality and the behavior of the 
derivatives of the diffusion coefficients changes drastically 
compared to the anisotropic cases. Looking at Eq. 
and taking into account that rotational symmetry implies 
Gxzx = Gyzy, we obtain A^^ = Aj, = -Gz - G^zx/'^- At 
variance with the behavior of the derivatives of the three 
energy gaps (compare Fig|21 and Fig^, all the deriva- 
tives of the three accelerations A^ — Ay and A^ diverge 
when approaching the critical point with a behavior anal- 
ogous to that of the derivative of the energy gap AE^ (see 
Eq. 0): 



dAa 

~dX 



d'^{x'-xir"' 



(7) 



with d'^^d'y^ -l/\/27r2 and 4 = -\/2/7r. 

In conclusion, we have introduced and studied ob- 
servable quantities, either static or dynamic, associated 
with single qubit unitary operations, that are able to 
characterize quantum phase transitions in spin systems. 
This provides an alternative, quantum-informatic point 
of view complementary to the traditional one in terms 
of correlation functions. Moreover, at variance with the 
characterization of quantum critical points based on mea- 
sures of entanglement, the quantities that we have intro- 
duced are, in principle, amenable to direct experimental 



observation. It is matter for future studies to establish 
whether this formalism can be useful for the character- 
ization of other quantum critical systems and models, 
and in different research contexts. In this last respect, it 
would be particularly interesting to investigate possible 
connections with quantum information transfer, informa- 
tion speed, and quantum speed limits in critical systems 



tion speed. 
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